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Abstract - Recent confocal experiments on colloidal solids motivate a fuller study of the projec- 
tion of three-dimensional fluctuations onto a two-dimensional confocal slice. We show that the 
effective theory of a projected crystal displays several exceptional features, such as non-standard 
exponents in the dispersion relations. We provide analytic expressions for the effective two- 
dimensional elastic properties which allow one to work back from sliced experimental observations 
to three-dimensional elastic constants. 



Optical techniques, including scattering and micros- 
copy, have long been used to extract detailed static and 
dynamic information from soft condensed matter systems. 
They are in many ways complementary - scattering being 
most suitable for examining the fluctuations in Fourier 
space El, giving information on the mode structure for 
uniform systems; microscopy gives the best account of the 
real space structure of a medium [2j and is particularly use- 
ful in the study of heterogeneous material properties. Re- 
cently several experimental groups have studied the prop- 
erties of colloidal crystals using video and confocal mi- 
croscopy at interfaces [3jj5] and in full three-dimensional 
samples |6j. Using computer analysis one combines the 
advantages of scattering and direct observation: One can 
observe a carefully chosen part of an experimental system 
and then study the mode structure of thermally excited 
fluctuations. It is particularly interesting to work deeply 
within a three dimensional sample [7 - 10 because one can 
be sure that surface perturbations to the properties are 
weak. 

However, the measurement of fluctuations rather than 
just the mean positions of the particles is technically diffi- 
cult: The scan must be fast in order to freeze the particle 
motion during the acquisition of a frame. Because of the 
constraints several groups have chosen to perform obser- 
vations on single confocal slices, rather than scanning the 
full three dimensional volume. The microscope resolves 
the motion of the colloidal particles within a single plane 
of the crystal structure-typically the plane (1,1,1). The 
sample is filmed for several minutes and the matrix of cor- 



relations is generated within the slice. Many thousands of 
frames are required in order to generate good fluctuation 
statistics. Experimentalists have studied weakly [10] or 
strongly disordered [8],[TT] materials with the hope of bet- 
ter understanding glassy dynamics, and characterizing the 
disorder via the projected fluctuation properties including 
the spectrum [8] , the eigenvectors [7] and the effective dis- 
persion relations [lO] . We show here that even in the case 
of ordered elastic materials a number of interesting fea- 
tures appear. Exceptional behaviour in disordered mate- 
rials should thus be defined with respect to the conclusions 
we present here. Observation of non-standard exponents 
for correlations within the slice can not be interpreted im- 
mediately as evidence of exotic or glassy behaviour. 

Our paper aims to calculate the effective theory which 
best describes the fluctuations of such a two-dimensional 
slice of a larger three-dimensional sample, in order to 
be able to easily work back from the observed two- 
dimensional correlations to three-dimensional material 
constants. The link between these elastic constants is 
given in Eqs. ( 14 1 below, which is the main result of the 
paper. The calculation of the projected properties will be 
based on the following principles: At large length scales 
a colloidal crystal is described by an effective elastic the- 
ory [i~||4{[T2] ; such an elastic theory leads to Gaussian fluc- 
tuations. Gaussian systems are rather special in that they 
allow one to exactly trace out degrees of freedom leading 
to a new effective theory which is also Gaussian in nature. 
This effective theory requires, however, effective renormal- 
ized couplings. Arguments due to Peierls [13] immediately 
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show that the final theoretical description of an elastic slice 
must be unusual: A two-dimensional elastic solid has di- 
verging fluctuations in the positions of individual particles, 
whereas the slice of the three-dimensional solid must have 
bounded fluctuations. The effective theory describing the 
slice can not be a simple variation on standard elasticity 
theory, we understand at once that anomalous dispersion 
relations are to be expected. 

Indeed we will show that while standard elasticity gives 
rise to a scaling of the elastic energy in q 2 as a function 
of the wavevector, q, the projected effective theory for the 
confocal slice is characterized by an effective dispersion 
relation in |q|. We give explicit analytic expressions for 
the prefactors in the dispersion law as a function of the 
three-dimensional elastic moduli, see Eqs. (14 1 below, al- 



lowing one to deduce the three-dimensional properties from 
measurements on the slice. The procedure described in 
this letter, albeit applied to elasticity, is sufficiently gen- 
eral that similar scalings in |q| should be observed in very 
different projected systems, such as in Stokesian hydrody- 
namics. 

We begin by linking the fluctuations in an elastic solid 
to the static elastic Green function of the medium. We 
then show how this Green function can be projected into a 
single layer, to produce an effective theory for the observed 
slice. We have performed extensive numerical simulations, 
which we compare with the analytic theory. 

Let us now consider a three-dimensional cubic crystal. 
Under small deformations the system is characterized by 
the displacement vector Ui and the symmetric homoge- 
neous tensor of displacement gradients My 14 . The elas- 
tic energy is then written as a quadratic form in mj which 
respects the cubic symmetry of the crystal. This quadratic 
form is related to the elastic matrix 
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Afe(k) = XS^Sm + fJ-(Si k dji + 5uS jk ) + vS ijk i kjk h (1) 

with Lame constants A, /i and anisotropy f FJ The hat 
denotes a Fourier transform and summation over repeated 
indices is assumed throughout the paper. The tensor S = 
ELi eWeWeWeW, with unit vectors parallel to the 
cubic axes of the crystal. The Green function of the static 
clastic problem is then the inverse of the elastic tensor, 



A,(k)Mk)=<5, 



(2) 



One expresses the free energy in terms of the displacement 
field 

F[fi] = i^fi i (k)6y(k)fi i (-k). (3) 
k 

If the crystal is studied at a finite inverse temperature /3, 
this implies that the correlation in the fluctuation ampli- 



1 Note that if the reference configuration of the crystal is under 
external stress, this stress appears explicitly in the elastic tensor [15] . 
A hard-sphere crystal is always under external pressure to be me- 
chanically stable. The elastic constants in Eq. |TJ implicitly contain 
this pressure correction. 
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tudes is given by 
(u,-(k)%(-k)) = 

^ ./1R ■ 

(4) 

For each wavevector k, D is a 3 x 3 matrix with eigenvalues 
c?i(k) where the subscript i indicates a polarization state. 
Following a convention usual in the experimental litera- 
ture [i][7] , we define the auxiliary variable uf (k) = di (k) FJ 

Instead of the full crystal, we now consider a crystal 
layer observed in a confocal microscope. In the following, 
Q a /3 is the Green function reduced to two dimensions, 
a, (3 £ {1,2}, and x,q £ R 2 are direct and reciprocal 
vectors in reduced space, whereas their three-dimensional 
counterparts are denoted r,k£ R 3 . Of course, neglect- 
ing the third dimension does not change the correlations 
within the layer; the real-space Green functions of the pro- 
jected and of the full problem are the same. 

From the Green function in the reduced space we then 
perform an inverse, two-dimensional, transform to find the 
effective dispersion relation for the observed slice. We wish 
to describe the fluctuations in the two-dimensional plane 
using a closed theory, calculating the two-dimensional 
equivalent of the matrix D. The calculational route that 
we will follow is 



— > Qofl(x) A Q af3 (q) — > D a0 (q), (5) 

where Ti is a ^-dimensional Fourier transform. 

The Green function C?jj(k) in three-dimensional recip- 
rocal space follows from the inversion of the elastic ma- 
trix Dij in Eq. Q. The result will have the following 
tensorial form, 



&i(k) 
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A" } ^f+A>"A ij (k) 



(6) 



where the prefactors A' , A" , A'" are cubic scalars and are 
not required to be isotropic; they may depend on the ori- 
entation of k. The first two terms in the brackets are the 
nearly isotropic part, while the third term Aij comprises 
further anisotropic properties. An explicit form of the 
latter can be obtained either by direct inversion of Dij us- 
ing the Sherman-Morrison formula or in a more elaborate 
way using bases for the space of all cubic tensors. The 
scalar prefactors are determined from the linear system 
of three equations which are obtained after multiplication 
of Eq. ([2]) with Sij, kikj/k 2 , and Aij, respectively. The 
Fourier transform has a similar tensorial structure, 
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(7) 



Again, the cubic scalars B',B",B"' depend on the orien- 
tation of k and are obtained as the solution of a linear 
system of equations. 



"Notice that ui is not the angular velocity of a wave. 
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The reduction to the two-dimensional Green func- 
tion Q aj a(x) can now be performed in real space. We 
firstly recognize that the result will be nearly isotropic, 
since the crystal plane we project on is a hexagonal lat- 
tice 114] . The Green function comprises two tensorial 
part^P 



Qa/3(X) 
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C'(x)A Q/J + C"(x)' 



(8) 



where again the hexagonal scalar prefactors may depend 
on the orientation of x. They are determined by pro- 
jection of the three-dimensional Green function onto the 
two-dimensional subspace. For doing this, we choose an 
orthonormal basis (rW , r^ 2 \ N), aligned such that N = 
(1, 1, l)/\/3 is orthogonal to the plane we project on. We 
now interpret x as a three-dimensional vector, denoted by 
an overbar, x = Y^ a =i x a Y^ a \ The two-dimensional iden- 
tity tensor S a /3 then becomes [Sij — NiNj] in three dimen- 
sions. Reduction of Eqs. (|8| and Q by 5 a p and x a xp/x 2 
and by their three-dimensional counterparts, allows to de- 
termine the prefactors in Eq. ([8]), 



2C"(x) + C"(x) 
Attx 



C"(x) + C"(x) XiX jn 
4^ = ^tM*) 



= [5 ij -N i N j ]g ij {yi), (9a) 

(9b) 



In the same way as in three dimensions, we obtain the re- 
duced Green function in reciprocal space using the ansatz 



Q«/3(q) 



i 

2q~ 



£'(q)<5 a/3 +.E"(q) 



q a qp 



(10) 



During the manipulations from Eq. ([6]) to Eq. ( 10 ) 
the scalar prefactors inherited the nontrivial vector- 
dependence from each other. It is important to notice 
that this dependence is restricted to the orientation of 
the vectors, since we regard only the long- wavelength 
limit, in which the elastic moduli in Eq. ([T]) are constants. 
The dependence on the norm of the vector is written 
out explicitly in Eqs. (J6J)— ( 10 (. In particular, the two- 
dimensional Fourier transform of Q(x) ~ 1/x led to the 
scaling Q(q) ~ l/q. 

Using only symmetry arguments, the scalar prefactors 
in Eq. ( |To| cannot be expected to be fully isotropic or 
even constants. In the numerical simulation described 
below, we observe however that in the limit of small q 
their angular dependence is negligible. Unfortunately, the 
direct-space cubic Green function Gij(v) cannot be calcu- 
lated explicitly, except for a few directions of higher sym- 
metry [16] . A more practical way is to approximate the 
cubic Green function by an appropriate isotropic one. Fe- 



dorov 17 18 provided an optimal way to do this, based 
on slowness curves. He proposed the effective isotropic 
moduli 

A = AH — , and fx = fi-\ — , (11) 
5 5 



3 If the external stress mentioned in note [T] were anisotropic, a 
third tensorial term would be allowed by symmetry. 



which give the optimal three-dimensional Green function 
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A 



fj, 



/2(A + 2/z) 



(12) 



The scalar prefactors in this and in the other Green func- 
tions are constants; A' and A" in terms of /2, A can be 
read off Eq. (pj); A'" = B'" = by isotropy; for the oth- 
ers, we find B' = C = A' + A" /2, B" = C" = -A" /2, 
E' = A', and E" = A"/2. The latter two give the isotropic 



projected Green function from Eq. (10 1, which we choose 



here to write in its longitudinal/transverse form 



1 A + 3/2 q a qp 
9 4/2(A+2/t) q 2 



(13) 



The effective elastic matrix D a p(q) of the projected two- 
dimensional slice, which is the inverse of this Green func- 
tion, thus predicts the following effective dispersion rela- 
tions: 



2/2 q (transverse) 



2 4/2(A+2/2) 
Wii = — = q 



A + 3/2 



(longitudinal). 



(14) 



These equations explain the anomalous (non-Debye) scal- 
ing of the projected fluctuations, where the auxiliary vari- 
ables u) are proportional to ^fq at long wavelengths, rather 
than q in standard elastic theory. The scaling is the same 
for all branches. The prefactors are here given explicitly in 
terms of isotropic approximations of the cubic elastic mod- 
uli for the hexagonal (1,1,1) plane in the face centered cu- 
bic crystal. For other modes and other orientations of the 
crystal explicit expressions are difficult to obtain, but nu- 
merical methods can help 19 . These simple expressions 



can be used by the experimentalist to obtain quantitative 
information on the three-dimensional elastic properties by 
studying only the two-dimensional slice. 

The distinction between the longitudinal and transverse 
branches has already been noted in the experimental re- 
ports [9][l0], which however do not try to fit the experi- 
mental prefactors in order to measure three-dimensional 
properties. 

To better understand the limits of our theoretical cal- 
culation we performed a molecular dynamics simulation 
with N = 4, 147, 200 particles organized in a face-centred 
cubic crystal, with a volume fraction </> = 0.57. We used 
event driven methods 
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because of their efficiency and 
also their long-time stability. The periodic simulation box 
contained TVi x iV 2 x iV 3 = 160 x 160 x 162 particles and 
had a skew shape, aligned with the Bravais lattice. This 
choice allows easy data analysis since the simulation box is 
also aligned with the slice direction (1, 1, 1), and it avoids 
ambiguities in the definition of the reciprocal vectors k. 

In order to study the mode structure of fluctuations we 
calculated and recorded the time average (ui(k)uj(— k)) in 
three dimensions and deduced the polarization eigenstates 
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Figure 1: Inset: Full, three-dimensional dispersion curves uJi/k 
used to extract the three elastic constants of a cubic crystal. 
We plot the longitudinal and transverse modes in the direc- 
tions (1,0,0) (red), (1,1,0) (blue) and (1,1,1) (black). The 
upper curves are longitudinal. For the orientations (1,0,0) 
and (1, 1, 1) the transverse modes are degenerate. Main curves: 
Dispersion relations in the confocal cut evaluated for the direc- 
tions (1,0), (2,1), (3,1), (4,1) with respect to the hexagonal 
Bravais lattice. The lines are the analytic prediction for the 
long-wavelength limit, Eqs. (14 1 and (111. 58,000 hours of 
computer time. 



by diagonalizing the three dimensional matrix of corre- 
lations measured for each k. For the three-dimensional 
modes we find the expected scaling di(k) ~ k 2 so that 
Ui ~ k. Thus when we plot u>i/k as a function of k (in- 
set of Fig. [IJ we can relate the small wavevector intersect 
to the three-dimensional isothermal elastic constants in 
Eq.Q. The transverse and longitudinal dispersion curves 
are then uniquely identified by their degeneracjQ and they 
contain sufficient information to extract the three indepen- 
dent elastic constants of a cubic crystal (2l|. We find the 
numerical values A, /x, v = 42.8, 51.8, —53.8 Units are cho- 
sen such that diameter, mass of all particles, and kT are 
unity. 

We now repeat the analysis for the sliced fluctua- 
tions of the crystal, (u a (— q)up(q)), using the hexagonal 
(1, 1, 1) planes to reproduce the experimental situation of 
(7|[l0, 22 , and extract the corresponding eigenvalues and 
auxiliary variables. x',(q), from the resulting 2x2 matri- 
ces. The result is plotted in the main figure of Fig. [TJ We 
see that two branches are important at long wavelengths, 
and that as predicted in our analytic theory the disper- 
sion relation for the modes are of the form ojf(q) ~ q. 
Anomalous projected dispersion curves have already been 
observed experimentally in both disordered and crystalline 
materials [l0l[22] . Using the effective elastic constants of 
Eqs. (fir} together with the dispersion law eq. (14), we 




calculate the prefactors to this law and plot the results as 
the lines in Figure [I] The theoretical and measured curves 
agree to within 3%. We interpret these discrepancies as 



see page 353 of Ref. [15] 



Figure 2: Inset: full dispersion curves for a two-dimensional 
hexagonal system. The upper curves are longitudinal, the 
lower ones are transverse; directions: red (1,0), blue (2,1), 
green (3,1), black (4,1). Main plot: Effective dispersion re- 
lation of a one-dimensional slice in the (1,0) direction, com- 
pared to the theoretical prediction (solid line). Simulation of 
400 x 400 particles, (f> = 0.85. 



being due to the anisotropy of the crystal. We note that a 
face-centred cubic crystal with nearest-neighbour interac- 
tions is predicted in linear elasticity to have A//i = 1 and 
v/fi = —1, see Eq. (12.7) of [23]. Our simulations find 
X/fj, = 1.21 and u/fi = -1.26. 

The use of the Fedorov form for the effective elastic 
constants is an uncontrolled approximation. To see to 
what degree the difference between theory and simula- 
tion is due to this approximation we also performed cal- 
culations and simulations of a two-dimensional ensem- 
ble of disks assembled in a hexagonal lattice projected 
onto a line. Since the two-dimensional elastic theory is 
isotropic |14| one can perform the projection without Fe- 
dorov's approximation within linear elasticity. We mea- 
sure again the two-dimensional dispersion curves in the 
inset of Fig. [2j and find a similar effective theory with 
uj 2 = 4q/i(A + 2//)/ (A + 3//) for the projection, which is 
the longitudinal branch of Eq. ( |14[ ) . The theoretical value 
for the coefficient is again plotted as a line. We find that 
there is no visible difference between the theory and the 
simulations for the projected system. 

To conclude, we have shown that if we wish to describe 
an elastic slice observed in a confocal microscope as an 
effective medium we must introduce an effective disper- 
sion relation in |q| which is very different from that which 
occurs in a normal two-dimensional medium with local 
interactions. Indeed in real-space one is obliged to con- 
sider that the system has long-ranged effective interac- 
tions. These interactions allow one to avoid Peierls' result 
implying that a two-dimensional system should not dis- 
play long-ranged order, due to the long-wavelength diver- 
gence of the expression for the mean squared amplitude of 
positional fluctuations: (u 2 ) ~ Jl/<? 2 d 2 q. This logarith- 
mically diverging result is replaced by a regular expression 
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due to the change in the dispersion law. 

We have shown that to good accuracy one is able to 
relate the three-dimensional elastic constants and two- 
dimensional elastic behaviour. Thus observations in two 
dimensions can be used to deduce estimates of the three- 
dimensional constants. It will be particularly interesting 
in the future to study how disorder and glassiness modify 
these effective properties j8j. It is interesting to note that 



24 where the 



such an energy function in |q| was found in 
spreading of a droplet was expressed as the effective dy- 
namics of a contact line. Again we are in the presence of 
a system projected to lower dimensions. 
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